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ABSTRACT 

The Kg band differential star count of the Two Micron All Sky Survey 
(2MASS) is used to derive the global structure parameters of the smooth com- 
ponents of the Milky Way. To avoid complication introduced by other fine struc- 
tures and significant extinction near and at the Galactic plane, we only consider 
Galactic latitude \b\ > 30° data. The star count data is fitted with a three- 
component model: double exponential thin disk and thick disk, and a power law 
decay oblate halo. Using maximum likelihood the best-fit local density of the 
thin disk is uq = 0.030 ± 0.002 stars/pc'^. The best-fit scale-height and length 
of the thin disk are Hzi = 360 ± 10 pc and Hri = 3.7 ± 1.0 kpc, and those 
of the thick disk are Hz2 = 1020 ± 30 pc and Hr2 = 5.0 ±1.0 kpc, the local 
thick-to-thin disk density ratio is /2 = 7± 1%. The best-fit axis ratio, power law 
index and local density ratio of the oblate halo are k = 0.55 ± 0.15, p = 2.6 ± 0.6 
and fh = 0.20 ± 0.10%, respectively. Moreover, we find some degeneracy among 
the key parameters (e.g., rio, Hzi, J2 and Hz'i)- Any pair of these parameters are 
anti-correlated to each other. The 2MASS data can be well-fitted by several pos- 
sible combinations of these parameters. This is probably the reason that there 
is a wide range of values for the structure parameters in literature similar to 
this study. Since only medium and high Galactic latitude data are analyzed, the 
fitting is insensitive to the scale-lengths of the disks. 

Subject headings: Galaxy: general - Galaxy: stellar content - Galaxy: structure - 
Galaxy: fundamental parameters - infrared: stars 
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Introduction 



In the eighteenth century, the famous astronomer Wilham He rschel showed u s the 



17851). The 



powerful method of star count to understand our own Milky Way f lHerschel 
technique has been used since then by generations of astronomers. With great improvement 
on data collection over the years, more and more details of our Milky Way were unfolded. 
However, the characteristic scales of smooth Galactic structures (i.e., disks and halo) 
obtained by previous studies does not converge to a common value as an outcome of 
improving data collection (see Tabled]). The spread of values is attributed to the degeneracy 
of Galactic model paramet ers (i.e., same star count data couM be fitted equally well by 



differ ent Galactic models) f lChen et al 



2001 



Siegel et al 



2002; 



Juric et al 



2008 



Bilir et al. 



20081 ). This is due to the diffe rent sky regions and limiting magnitudes (i.e., limiting 



volumes) used in these studies (ISiegel et al. 



Juric et al. 



20081 ). On the contrary. 



2002 



Bilir et al. 



Karaali et al. 



(120081 ): 



2004; 



3ilir et al. 



Yaz fc Karaalil (120101 ) did not show 



2006ab 



such degeneracy in determining the Galactic model parameters. This controversy is still 
actively debated. Therefore, systematic all sky surveys with deeper limi ting magnitude 
and w ider sky region, such as the Two Micr on All Sky Survey (2MASS; 



Skrutskie et al. 



20061 ). the Sloan Digital Sky Survey (SDSS; 



Yor. 



c et al. 



elescope fc Rapid Re sponse System (Pan-Starrs: 



2000 



, the Panoramic Survey 



Kaiser et al. 



( jPerryman et al. 



20021 ) and the GAIA mission 



200ll ) could provide a good opportunity for us to study our Galaxy from 



a global perspective. Free from limited sky fields, astronomers can acquire many more 
information from the stellar distribution of these surveys. 

On Galactic structu re study, besides the simple and smooth two-component model 



( JBahcall fc Soneira 



19801 ) or the three-component model (JGilmore fc Wyse 



more structures have been discovered, such as inner bars in the Galactic center (lAlves 



Hammersley et al. 



2000; 



van Loon et al. 



2003 



Nishiyama et al. 



2005 : 



1985|). many 



2000 



Cabrera-Lavers et al. 
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2006 



2008 ) , and flares and w arps ( iLopez-Corredoira et al.ll2002l : iRobin et al.ll2003t iMomany et al. 



Reyle et al 



20091 ) . which has been contributed t 



with Galactic longitude (ICabrera-Lavers et al 



2007 



le variation of d i sk model pararn eters 



Bilir et al 



20081 



20101). Moreover, the ove rdensities in the halo, such as Sagittarius 



Trian gulum- Andromeda (JRocha-Pinto et al 



2004 



Maiewski et al. 



20081). and in the ou t er disk, such as Canis M ajor (JMartin et al. 



( iNewberg et al. 



2002; 



Rocha-Pinto et al 



Yaz Sz Karaali 



(IMaiewski et al. 



2004h. Virgo ( jjuric et al. 



2003h. 



20041 ) ■ Monoceros 



20031 ). show the complexity of the Milky Way. 



The formation history of our Galaxy is more complicated than what we thought previously. 



On stellar luminosity function study, a recent study by IChang et al.l (120 



of) us ing the 



Kg band star count of 2MASS point source catalog (2MASS PSC; ICutri et al.ll2003l ) verifled 
for the flrst time the universality hypothesis of the luminosity function (i.e., the common 
practice of assuming one luminosity function for the entire Milky Way). 

We are interested in the global smooth structure of our Galaxy. In view of the 
existing flne structures (e.g., flares, warps, overdensities. . .etc.), the structure of the smooth 
components can be determined either by including these fine structures in a grand Galaxy 
model or by avoiding the sky area "contaminated" with these features. The first method 
demands a complex model, which involves many mor e structure parameters, and ne eds 



Lopez- Corredoira et al. 



(120021) is a 



high computing power to accomplish the fitting task, 
good example of this method using 2MASS data. The second method is clearly simpler but 
needs justifications. We observe that (1) the fine structures (e.g., inner bars, fiares and 
warps), which have observable contribution on 2MASS star count data, are all confined in 
the Galactic plane region. (2) The overdensities or substructures in the outer disk region 
and the halo are difficult to identify in general. Their c ontrib ution is negligible in the 



2MASS star count data. Here we quote 



Majewski et al. 



(120041 ) on halo substructure: "This 



substructure is typically subtle and obscured by a substantial foreground veil of disk stars. 
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eliciting its presence requires strategies that optimize the substructure signal compared to 
the foreground noise." 

We prefer the second method in this paper and only use Galactic latitude \b\ > 30° 
2MASS Kg band star count data to obtain the structure parameters of a three-component 
model. In addition, the influence of the near infrared extinction in these regions is small. 
This also allows us to use a simpler extinction model for correction. We describe our model 
in section 2 and the analysis method in section 3. Section 4 provides the results and a 
discussion. 

2. The Milky Way Model 

We adopt a three-component model for the smooth stellar distribution of the Milky 



Way. It comprises a t hin disk, a thick disk and an oblate halo ( jBahcall fc Soneira 



1980 



Gilmore fc Reid 



19831). The total stellar density n{R, Z) at a location [R, Z) is the sum of 



the thin disk Di, the think disk D2 and the halo H, 



n{R, Z) = no [D,{R, Z) + D^iR, Z) + H{R, Z)\ , (1) 

where R is the galactocentric distance on the Galactic plane, Z is the distance from the 
the Galactic mid-plane and no is the local stellar density of the thin disk at the solar 
neighborhood. 

The stellar distribution of the thin disk D\ and the thick disk D2 decreases exponentially 
along R and Z (the so called double exponential disk), 

{R-R^) {\Z\-\Z^\) 



Di{R,Z) = fiQxp 



(2) 



Hri Hzi 

where {Rq,Zq) is the location of the Sun, !!„ is the scale-length, H^i is the scale-height, 
and fi is the density ratio to the thin disk at the solar neighborhood. The subscript z = 1 
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stands for t he thin dis k (thus /i = 1) and z = 2 for thick disk. We adopted Rq = 8 kpc in 



our model (JReid 



19931) 



The halo is a power law decay oblate spheroid flattening in the Z direction, 



H{R,Z) = h 



R^ + (Z/k) 



-p/2 



(3) 



Rl + {z^/^y 

where k is the axis ratio, p is the power index and fh is the local halo-to-thin disk density 
ratio. 



Chang et al.l (120 lOf ) showed that the whole sky Kg band luminosity function can be 
well approximated by a single power law with a power law index 7 = 1.85 ± 0.035, a 
bright cutoff at Mb = -7.86 ± 0.60 and a faint cutoff at Mj = 6.88 ± 0.66. We adopt this 
luminosity function in the following analysis. In normalized form, it is 

2 log, 10 (7-1) 



^(M^ 



;^q2(7-1)M/5 _ 



(4) 



^ 5 ri02(7-l)M^/5 _ ]_02(7-l)Mj5l 

Note that ip{M) includes all luminosity classes. In our analysis the observing magnitude 
range is 5 < -ft's < 14 mag. The corresponding distances of bright cutoff, faint cutoff and 
Mk^ = are 3.4 kpc to 216 kpc, 4 pc to 265 pc and 0.1 kpc to 6.3 kpc, respectively. 

Although we only use NIR data in the medium and high Galactic latitude regions, we 
still need to corre ct possible interstellar extinction. We adopt the new COBE/IRAS result 



Chen et al 



(ISchlegel et al. 



1999h a nd convert it to the Kg band extinction by AxjE{B — V") = 0.367 



19981 ). This extinction model is then applied to our simulation data. The 



extinction values of most of our analyzed regions are A^^ < 0.03. 
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3. The Data and Analysis Method 



3.1. The 2MASS Data 



The 2MASS Point Source Catalog (2MASS PSC, ICutri et al.ll2003h is employed to 
carry out the Kg band differential star count for the entire Milky Way. We div ide the whole 
sky i nto 8192 nodes according to level 5 Hierarchical Triangular Mesh (HTM 



Kunszt et al. 



200ll ). The level 5 HTM samples the whole sky in roughly equal area with an average 
angular distance about 2 degrees between any two neighboring nodes. The amount of 
stars within 1 degree radius of each node (i.e., each node covers vr squ are degree) i s the n 



retrieved with a bin size -ft's=0.5 mag via 2MASS online data service (jCutri et al 



2003 . 



catalog). Our selection criterion is: the object must be detected in all J,H,Ks bands and 
has signal-to-noise ratio > 5. Since the limiting magnitude of 2MASS Kg band is 14.3 ma g, 
which has 10 signal-to-noise ratio and 99% completeness (see Table 1 in 



Skrutskie et al. 



20061 ). and Kg < 5 mag objects have relatively large photometric error, we only compare 
2MASS data with our simulation data from i^g = 5 to 14 mag. 

In order to minimize the effects coming from the close-to-Galactic-plane fine structures 
(e.g., flares, warps, arms, budge and bars. ..etc., which have considerable contribution 
to the star count data), and the relatively complex extinction correction at the low 
Galactic latitude region, we avoid low galactic latitudes, and only consider data in 
Galactic latitude \b\ > 30°. Although several overdensities in the halo (such as Sagittarius, 
Triangulum-Andromeda, Virgo. . . etc.) were identified, they cannot be picked up from the 
overwhelming foreground field sta rs on star count data w ithout additional information (e.g.. 



color, distance, metallicity. . . etc.. 



Majewski et al. 



20041 ). Therefore, their contribution to 



the 2MASS Kg band differential star count is negligible and will not affect our result. We 
also exclude the areas around Large and Small Magellanic Clouds for their significant stellar 
population. 



3.2. The Analysis Method 

The Maximum Likelihood Method (Bienayme et al. 1987) is applied to compare the 
Kg band 2MASS differential star counts and the simulation data to search for the best-fit 
structure parameters of the three-component Milky Way model. Our fitting strategy is as 
follows: 

1. Take Rq = 8 kpc; 

2. choose one Zq and work out the maximum likelihood value by fitting the 9 parameters 
(no, Hzi, Hri, /2, Hz2, Hr2, fh, K,p); 

3. repeat step 2 for other Zq; 

4. pick the Zq corresponds to the maximum of the maximum likelihood values in step 3; 

5. repeat steps 2 to 4 for finer grid size of the 9-parameter fit and a narrower range of 
Zq around the one found in step 4; 

6. the uncertainty is estimated by adding Poisson noise on the simulation data to see 
how the likelihood varies. The difference of the likelihoods of 500 realizations of the 
same model, differed by the Poisson statistics only, gives a range of likelihood around 
the maximum likelihood that defines the confidence level. 

Table [2] lists our searching parameter space and the finest grid size we used. The key 
parameters in our study are uq, H^i, /2 and Hz2 (see Eqs. (l)-(2)). The first two play a 
primary role on the variation of the likelihood value and the latter two play a secondary 
role. The other five parameters Hri, Hr2, fh, i^ and p (see Eqs. (l)-(3)) are non-key 
parameters, which play a minor role and do not affect the likelihood value as much as the 
key parameters. 
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The Results and Discussion 



Table [3] lists our best-fit results with the corresponding uncertainties. Fig. [T] shows 
contour plots of likelihood against different pairs of parameters. The contour changes 
dramatically along the key parameters uq, Hzi, /2 and Hz2, but relatively mild along other 
parameters Hri, Hr2, fh, P and k. This indicates the importance of the key parameters in 
determining the best-fit result. 

Degeneracies exist between some pairs of key parameters, such as (no,Hzi), (no,/2), 
{Hzi, /2). . .etc. Here degeneracy means that the likelihood value stays almost the same 
when the pairs of parameters change together in a particular way. Similar degeneracy 
between the local th ick-to-thin d i sk ra ti o f-?, and the scale- h eight of the thick disk ff^^ 



has be en reported in 



Chen et al. 



fl200lh 



Siegel et al. 



(120021 ): 



Juric et al. 



mm; 



Bilir et al. 



(120081 ). Consequently, it is possible that different combinations of parameters can be 
regarded as 'acceptable' fitting. For example, if we choose a higher local stellar density 
uq, then we can pick a smaller H^i such that the likelihood value is very close to the 
maximum likelihood value and assign it as the 'best-match' scale-height of the thin disk. 
Therefore, a thin light color diagonal strip shows on the Uq against H^i contour plot 
in Fig. [U Besides, similar trends happen in other pairs of key parameters. When one 
parameter of the pair is higher, we can get a similar likelihood value by lowering the 
other parameter (see the corresponding two-parameter contour plots of key parameters 
in Fig. [T]). This anti-correlation is not unexpected. The number of stars along the line 
of sight in the model increases when any one of the key parameters increases. Thus for 
a given observed number of stars, an increase in one key parameter can be compensated 
by a decrease in the other. Perhaps this is the reason that our best-fit scale-height of 
the th in disk, H^^ = 360 pc, i s sorn ewhat larger than the reported values, H^i = 285 



pc, m 



Lopez-Corredoira et al. 



(I2OO2I ) study (they also use star count of 2MASS to obtain 
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Galactic model parameters). Our best-fit local stellar density from Kg =-8 to 6.5 mag is 
riQ ~ 0. 030 star/pc'^, which is about half of ~ 0.056 star/pc'^ of the corresponding value 



cited in lEaton et al 



value cited in 



m 



19841). and a bit 



Lopez-Corredoira et al. 



ower than ~ 0.032 star/pc^ of the corresponding 



( I2OO2I ). As a result, our fitting tends to choose a 



larger scale- height of the thin disk. I 



with that of 



Lopez-Corredoira et al. 



we force our local stellar density to be comparable 



(J2OO2I ). then the corresponding 'best-fit' scale-height of 
the thin disk would be ~320 pc, which is closer to their result. If we choose even higher 
local stellar density, then the 'best-fit' scale-height of the thin disk would be made between 
200 to 300 pc, which is similar to the most of recent studies (see Tabled]). Moreover, we 



do not apply binarism correction in our analysis, and it has been show n that scale- 



and scale-height might be underestimated without bina rism correction (jSiegel et al. 



Juric et al. 



2008 



Ivezic et al. 



2008 



Yaz&Karaali 



ength 



2002 



2010h . 



Girardi et al. 



Wainscoat et al. 



20051) o r luminosity 



1992) do exist, but 



For our purpose, we deem that in order to lift the degeneracy it is crucial to have a 
reliable near infrared luminosity function by observation or a near infrared local stellar 
density. Unfortunately, a systematic study in this direc tion is yet to come . Some related 
studies, such as synthetic luminosity function (see e.g., 
function transformed from optical observation (see e.g., 
some uncertainties still need to be settled (e.g., the initial mass function, mass-luminosity 
relation for NIR and color transformation between different wavelengths. . .etc.). Once the 
'true' local stellar density is known, the 'true' structure of the Milky Way would be revealed. 

In order to see how good the agreement between 2MASS data and our best-fit model, 
we show an all sky map of the ratios of observed to predicted integrated star count from 
Ks =5 to 14 mag for each node as a function of position on the sky in Fig. [21 and the 
color indicates the values of the ratios. We do not see obvious deviation in the Galactic 
latitude |6| > 30° areas, but only in the Large Magellanic Cloud and the Small Magellanic 
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Cloud areas. For comfirmation, we plot integrated star count from Ks =5 to 14 mag 
along Galactic longitude and latitude for \h\ > 30° areas in Figs. El&lll We see 2MASS 
data and our best-fit model agree well and only some small deviations in the nodes at 
the anti-Galactic center h ~ 30° areas. The significant spikes in Figs. [3] & H] are due to 
the populations of the Large Magellanic Cloud and the Small Magellanic Cloud. For 
testing how these small deviations affect the key parameters selection, we exclude these 
small deviations and re-analyze the data with fixed non-key parameters. It shows only 
the scale-height of thick disk shifts slightly but still within our error estimation. Besides, 
we do not see any significant differences of reported overdensities in sky regions with 
\h\ > 30°. Thus, we conclude that the three-component model can describe the Milky 



Way structure sufficie ntly well for 



luminosity function of 



l igh G alactic latitude regions and the single power law 



Chang et al. 



purpose is to search for the global Galactic mode 
best-f it result for each node individually as what 



(120101 ) is a good approximation as well. Since our main 
parameters, we do not try t o explore th e 



Cabrera-Lavers et al 



(120071 ) 



Bilir et al. 



( I2OO8I ): lYaz fc Karaalil ( I2OIOI ) have done in their studies to seek the variations of disk model 
parameters with Galactic longitude. Instead, we treat the differences of 2MASS data to 
our best-fit model as deviations from a global smooth distribution. We believe that similar 
variations in disk model parameters will be obtained if we take similar analysis procedure 
(i.e., searching best-fit parameters for each node), but this is beyond the scope of this work. 
Because we do not consider flares, warps and other overdensities in our model, there are 
some discrepancies between 2MASS data and the model in the low Galactic latitude regions 
(see Fig. |2]). These fine structures make the star distribution more fluffy in the vertical 
direction toward the edge of the Milky Way. Hence the discrepancy between 2MASS data 
and the model increases vertically towards the anti-Galactic center region. In addition, the 
absence of Galactic bulge in our model contributes to the large discrepancy in Galactic 



center areas. The difference in the low Galactic latitude region needs more delicate analysis 
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to rectify (see, e.g. 



Lopez- Corredoira et al 



2002 



Moinany et al 



20061 ) 



4.1. Summary 

In summary, we set forth to study the global smooth structure of the Milky Way by a 
three-component stellar distribution model which comprises two double exponential disks 
(one thin and one thick) and an oblate halo. The Kg band 2MASS star count is used 
to determine the structure parameters. To avoid the complication introduced by the fine 
structures and complex extinction correction close to Galactic plane, we use only Galactic 
latitude \h\ > 30° data. There are 10 parameters in the model, but only four of them play 
the dominant role in the fitting process. They are the local stellar density of the thin disk 
77-0, the local density ratio of thick-to-thin disk /2, and the scale-height of the thin and 
thick disks H^i and Hz2- The best-fit result is listed in Table [31 In short the scale- height 
of the thin and thick disks are 360 ± 10 pc and 1020 ± 30 pc, respectively; the scale-length 
of the thin disk is 3.7 ± 1.0 kpc and that of thick disk is 5.0 ± 1.0 kpc (the uncertainty 
in scale- length is large because it is not very sensitive to high latitude data.) The local 
stellar density ratio of thick-to-thin disk and halo-to-thin disk are 7 ± 1% and 0.20 ± 0.10%, 
respectively. The local stellar density of the thin disk is 0.030 ± 0.002 stars/pc^. 

An all sky comparison of the 2MASS data to our best-fit model is shown in Fig. [2l 
A good agreement in the Galactic latitude |6| > 30° areas is expected from our fitting 
procedure. In low Galactic latitude regions, fine structures (such as flares, warps. . . etc.) 
increase the effective scale-height towards the edge of the Milky way. This is reflected in 
the fan-like increase in discrepancy towards the anti-Galactic center regions. 

Degeneracy (i.e., different combinations of parameters give similar likelihood values) is 
found in pairs of key parameters (see Fig. [1]). Thus different combinations of parameters 
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may fit the data almost as good as the best-fit one, and these are all legitimate 'acceptable' 
fitting in view of the uncertainty. Therefore, accompanying our lower local stellar density 
0.030 stars/pc'^ (from Kg = —8 to 6.5 mag) is a higher thin disk scale-height 360 pc. In the 
context of NIR star count, the NIR luminosity function or the NIR local stellar density is 
imperative to determine the scale-height and other Milky Way structure parameters. We 
hope that systematic study on the luminosity function and the local stellar density in near 
infrared will be available in the near future. 

We acknowledge the use of the Two Micron All Sky Survey Point Source Catalog 
(2MASS PSC). We would like to thank the anonymous referee, whose advice greatly 
improves the paper. CMK is grateful to S. Kwok and K.S. Cheng for their hospitality 
during his stay at the Department of Physics, Faculty of Science, University of Hong Kong. 
This work is supported in part by the National Science Council of Taiwan under the grants 
NSC-98-2923-M-008-001-MY3 and NSC-99-2112-M-008-015-MY3. 
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Table 1. Previous Galactic models. The parentheses are the corrected values for binarism. 
The asterisk denotes the power-law index replacing Re. References indicating with the 
original result table mean Galactic longitude or limiting magnitude dependent Galactic 

model parameters. 



H^l (pc) 


H,l 


(kpc) 


/2 (%) 


H,2 (kpc) 


H^2 


(kpc) 


fh(%) 


Rc(S) (kpc) 


H, 


Reference 


310 - 325 


3.8 

3.2 
2.5 
2.3 

4 

2.25 
2.8 
2 - 2 
1.97 
3.5 

2.1 
1.9 

1.65- 

2.15( 
1.00- 
3.7 


2.52 

;2.6) 
1.68 


0.0125-0.025 

0.02 

0.02 

0.0028 

0.016 

0.02 

0.015 

0.0224 

0.041 

0.019 

0.020-0.025 

0.019 

0.056 

0.074 

0.059 

0.061 

0.02 

0.065 - 0.13 

3.5 

0.06 - 0.10 

0.02-0.03 

0.07 

0.052-0.098 

0.11 

0.04-0.10 

0.033-0.076 

0.087/0.055 

0.01-0.07 

0.16/0.10 

0.062-0.145 

0.027-0.099 

0.055-0.151 

0.13(0.12) 

0.0640-0.0659 

0.07 


1.92 - 2.39 

1.45 

1.3 

1.9 

1.18 - 2.21 

1 

1.3 - 1.5 

0.95 

1 

0.9 

0.86 

1.6-1.4 

0.98 

0.76 

0.76 

0.91 

0.79 

1.26/1.29 

0.58 - 0.75 

0.86 

0.7 - 1.0 (0.9 - 1 

0.9 

0.64 

0.805-0.970 

1.06 

0.9 

0.84-0.87 

0.58/0.93 

0.6-1.1 

0.49/0.58 

0.80-1.16 

0.62-1.03 

0.55-0.72 

0.743(0.900) 

0.860-0.952 

1.02 


3.8 

4.3 
2.8 
3 
3 

3.5 
3.7 
.2) 3-4 

4.7 
3.04 

2.3-4.0 

3.261(3.600) 
2.65-5.49 


0.002 

0.0012 

0.0013 

0.001 

0.002 

0.001 

0.002 

0.0011 

0.0015 
0.0024 
0.0015 

0.0005 

0.0013 

0.0015 

0.002-0.003 

0.00125 

0.0002-0.0015 

0.0004-0.0006 

0.001 

0.00125 

0.0001-0.0022 

0.0007-0.0019 

0.0051 

0.0033-0.0039 

0.002 


3 

3.1* 

2.36 

2.9 

3 

2.7 

4 

2.67 

3.3 

2.44 - 2.75* 

2.69 

2.99* 

4.3 
3/2.5* 

2.77* 
2.6* 


0.85 

0.8 

0.85 

Flat 

0.9 

0.85 

0.84 

0.8 

0.48 

0.60 - 0.85 

0.84 

0.63 
0.55 

0.50 - 0.70 

0.5-0.6 

0.6 

0.6-0.8 

1/0.6 

0.06-0.08 
0.74 
>0.4 
0.45 

0.25-0.85 

0.53-0.76 

0.64 

0.0489-0.0654 

0.55 


Yoshii ^1982) 


300 


Gilmore & Reid (1983) 


325 
280 


Tritton & Morton (1984) 


125-475 

300 

285 


Robin & Creze a9861 
del Rio & Fenkart (1987) 


325 


Yoshii et al. (1987) 


249 
350 


Kuiiken & Gilmore a989'l 
Yamaeata & Yoshii (1992) 
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Table 2. The searching parameter space. 





Range 


Grid size 


Thin Disk 






Hrl 


1.0-6.0 kpc 


100 pc 


H,i 


200-450 pc 


10 pc 


no 


0.02-0.04 stars/pc^ 


0.002 stars/pc^ 


Zq 


11-35 pc 


2.5 pc 


Thick Disk 






Hr2 


3.0-7.0 kpc 


200 pc 


H,2 


400-1200 pc 


20 pc 


h 


0-20% 


1% 


Spheroid 






K 


0.1-1.0 


0.05 


P 


2.3-3.3 


0.1 


fh 


0-0.45% 


0.05% 



Table 3. The best-fit Milky Way model. 





Value 


Uncertainty 


Thin Disk 






Hrl 


3.7 kpc 


1.0 kpc 


H.l 


360 pc 


10 pc 


no 


0.030 stars/pc3 


0.002 stars/pc^ 


Zq 


25 pc 


5 pc 


Thick Disk 






Hr2 


5.0 kpc 


1.0 kpc 


H,2 


1020 pc 


30 pc 


h 


7% 


1% 


Spheroid 






K 


0.55 


0.15 


P 


2.6 


0.6 


fh 


0.20% 


0.10% 
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